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ABSTRACT 

A new, inexpensive polarimetric unit has been constructed for the Dominion 

Astrophysical Observatory (DAO) 1.8-ni Plaskett telescope. It is implemented 

as a plug-in module for the telescope's existing Cassegrain spectrograph, and 

p ■ enables medium resolution [R "^ 10, 000) circular spectropolarimetry of point 

sources. A dual-beam design together with fast switching of the wave plate at 

rates up to 100 Hz, and synchronized with charge shuffling on the CCD, is used to 

c^ ' significantly reduce instrumental effects and achieve high-precision spectropolari- 

metric measurements for a very low cost. The instrument is optimized to work 

^ ■ in the wavelength range 4700 — 5300 A to simultaneously detect polarization sig- 

00 ' nals in the H/3 line as well as nearby metallic lines. In this paper we describe the 

Li ■ technical details of the instrument, our observing strategy and data reduction 

O ■ techniques, and present tests of its scientific performance. 

rn 

pi '. Subject headings: Stars - Astronomical Instrumentation - Data Analysis and 

Techniques 

>: 
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C^ . 1. Introduction 

Spectropolarimetr y has been used for decades to detect magnetic fields in stars through- 



O 

H 



out the HR diagram (JDonati fc Landstreetl 120091 ). Our particular interest in this paper. 



however, is in instrumentation that enables the measurement of strong, large-scale magnetic 
fields in the chemically peculiar stars (Ap/Bp) of the upper main sequence. Some existing 
polarimeters permit high- resolution (defined here as i? > 40,000), high-precision spectropo- 
larimetry in metallic lines of such stars (for example, MuSiCoS and NARVAL on the 2 m 
Bernard Lyot Telescope at Pic du Midi, and ESPaDOnS on the 3.6 m Canada- France-Hawaii 
Telescope) but hydrogen lines, due to their large width, provide much lower accuracy. Low- 
resolution [R ^ 2,000) instruments such as FORSl on the 8 m Antu telescope of the VLT 
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( iBagnulo et al.l |2002| ) enable more precise spectropolarimetric measurements of hydrogen 
lines but at a cost of providing such low contrast for the metallic lines that they can not 
be fully utilized. Because of this, high-resolution spectropolarimeters are generally used to 
measure magnetic fields in Ap and Bp stars using only metallic (and perhaps helium) lines, 
while only hydrogen lines are measured with low-resolution instruments. 

However, for many years it has been known that a number of stars show quite significant 
discrepancies between magnetic field strengths derived from hydroge n lines versus values 



found using metallic lines, but with differ ent instruments (for example. iBorra fc Landstreet 



19801 : iBychkov. Bychkova. fc Madejl 119771 ). The differences can be from a few hundred G to 
a few kG, or more than 10% of the field strength since longitudinal magnetic fields in these 
stars can range from several hundred G to several kG. As an example, in 52 Her the average 
longitudinal field derived from metallic lines is 1500 G while the average field derived from 



hydrogen lines is 500 G (IBychkov. Bychkova. fc Madejl 119771 )! Unfortunately, it is unclear 



what part of the observed differences might be due to details of the instrumentation and 
data analysis techniques and what part may be produced by the very common and often 
very pronounced non-uniform horizontal and vertical abundance distributions of metals in 
the photospheres of the Ap and Bp stars. 

In an attempt to resolve this long-standing question we have investigated the possibil- 
ity of constructing a low-cost, medium- resolution spectropolarimeter that provides circular 
polarization measurements of hydrogen and metallic lines of Ap and Bp stars and hence 
enables the measurement of longitudinal magnetic fields with comparable precision using 
both diagnostics simultaneously. The low cost was dictated by the limited operations and 
development budget for the DAO telescopes. Despite this constraint, we set out to build 
an instrument that would let us measure stellar magnetic fields in relatively bright stars 
with typical accuracies of about 100 G and better than 50 G for the brightest targets. We 
will demonstrate later that such precisions are sufficient to routinely detect discrepancies 
between hydrogen and metallic line longitudinal field measurements and also show for the 
first time that for at least one object these discrepancies do not have an instrumental origin. 

A second science goal of the new polarimeter was to provide us with the ability to carry 
out (in a reasonable length of time) a large, sensitive search for globally-ordered magnetic 
fields in relatively faint and possibly rapidly rotating upper main sequence peculiar stars 
observable from the DAO. For hundreds of such poorly studied 6th and 7th mag Ap and 
Bp stars we therefore needed to obtain field measurements with an accuracy on the order 
of 200 G from observations acquired in two hours or less. Such measurements, obtained 
many times over the rotation period of the star, are sufficient to confidently detect and 
characterize the large-scale magnetic fields observed in many upper main sequence peculiar 
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magnetic stars. 

In the remainder of this paper we describe the optical and mechanical design of the 
spectropolarimeter 'dimaPol' that we have designed and constructed for use on the DAO's 
1.8-m Plaskett telescope for a very low cost of only several thousand dollars (including spare 
components and parts used only for testing). The polarimeter module can be switched at 
a rate of up to 100 Hz and is implemented as a plug-in module for the existing Cassegrain 
spectrograph on the DAO 1.8-m telescope. We describe our data reduction technique that 
not only enabled the use of low cost, off-the-shelf components in the polarimeter but also 
provides near real-time conversion of the spectropolarimetric data to a magnetic field mea- 
surement. A summary of tests we have carried out to investigate the optimal modulation 
rate of the polarimeter and the best observing strategy is provided, and we also present 
results of extensive new magnetic field measurements of well-known magnetic stars as well 
as non-magnetic standards and demonstrate how these new data compare very favourably 
to observations made with other generally much more costly spectropolarimeters. 



2. The DAO Liquid Crystal Spectropolarimeter 

The level of polarization observed in starlight is usually very weak. Numerous instru- 
mental effects can make these weak signals undetectable or lead to spurious detections. Most 
astronomical spectropolarimeters currently in operation use a dual-beam design to deal with 
this issue. The polarization of the input beam is transformed into an intensity difference 
between the two output beams. Since both spectra are taken under exactly the same con- 
ditions, the polarization signal is independent of sky transparency, seeing, or slit losses but 
not the instrumental systematics. A polarization modulator is therefore used in many spec- 
tropolarimeters to fight these possible instrumental effects by interchanging the two output 
beams so that they travel exactly the same optical path through the instrument. Usually 
the modulation is achieved by means of the mechanical rotation of an optical element but 
this can be a relatively slow process. Since at least two observations are now needed for each 
distinct configuration of the device and these are separated in time, guiding errors can now 
come into play. These can cause the intensity distribution of the input light beam across 
the entrance aperture of the spectrograph to change from one observation to the next which 
in turn can lead to a wavelength shift on the detector and hence may result in a spurious 
polarization signal. A change in brightness due to seeing and transparency variations and 
different gains in the two instrument c hannels can also generate a spurious signal (see the 



discussion of this effect in lPovel Il200ll ). This can be a limiting factor for the sensitivity of 



a spectropolarimeter unless the instrument is fed by an optical fiber (as in ESPaDOnS, for 
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Fig. 1. — A sketch of the optical layout of the polarimetric module. 



example) or some form of fast modulation is used. The spectrograph currently permanently 
installed on the DAO Plaskett 1.8-m telescope is fed directly from the telescope's Cassegrain 
focus. Because of this fact, as well as budget constraints, we confined our design consid- 
erations for a new DAO spectropolarimeter to those that made use of a fast modulation 
option. 

dimaPol is installed in the DAO 1.8-m telescope's //IS diverging beam immediately 
behind one of three possible circular entrance apertures to the Cassegrain spectrograph. The 
module consists of an achromatic quarter-wave plate, a ferro-electric liquid crystal (FLC) 
half- wave plate, a calcite beam displacer, and a mechanical shutter (Figured]). We discuss 
the individual optical components in detail below. The single-order grating spectrograph 
has an off-axis hyperbolic collim ator mirror, a plane diffraction grating, and a spherical 
camera mirror (JRichardsonlll968l ) . An 1800 lines mm~^ grating provides a spectral resolution 
of -R = 8, 000 to 15, 000 depending on the choice of entrance aperture. The 1752x532 SITe 
CCD with 15 /im pixels captures an approximately 260 A long spectral window centered on 
the H/3 line and the polarimeter optics have been optimized to work in this wavelength range. 
The spectral resolution is high enough to resolve spectral lines in slowly rotating stars but 
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at the same time is low enough so that any polarization signal in the broad hydrogen line, 
as well as broad metallic lines of rapidly rotating stars, can be precisely measured. 



2.1. Quarter- wave plate 

The quarter-wave plate in the polarimeter module converts left and right circularly 
polarized light in the input beam into linearly polarized light with orthogonal polarization 
directions. The achromatic quarter-wave plate in dimaPol was made by Bolder Vision Optik 
and is a multilayer polymer laminated between two optically flat windows. The retardance 
error is less than 0.004 A over the wavelength range where the polarimeter performance is 
optimized (Figure [2]). A polymer plate was selected because of its large acceptance angle, a 
requirement for use in the diverging input beam from the telescope. 
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Fig. 2. — Measured retardance of the quarter-wave plate used in dimaPol. The nominal A/4 
retardance is shown by a dashed line. Error bars are approximately the size of the data 
points. 
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2.2. Half- wave plate 



The quarter-wave plate is followed by a half-wave plate, a commercial FLC made by 
Displaytech Inc. (model LV2500-OEM). This is an electrically switchable retardation plate 
with a fixed retardation of approximately 180° or A/2 but the orientation of the optical 
axis can be switched between two distinct angles by changing the voltage applied across 
the liquid crystal. A commercial driver DR95 made by Displaytech Inc. is used to supply 
and switch the voltage. Switching commands are issued by the CCD controller. For FLC 
wave plates this switching can be very fast. In solar astronomy FLC wave plates have been 
used in a number of polarimeters, for example, the Zurich Im aging Polarimeter (ZIMPOL, 



Gisler et al.ll2004[ ). the SOLIS Vector- Spectromagnet ograph keller et all I2OO3]) LPS? & 



TIP at the solar telescopes on the Canary Isl ands (iMartinez Pillet et al.lll999l ). and the 



spectrograph of the Pic du Midi Turret Dome (iMalherbe et al.ll2004j ). Modulation rates of 



J 



up to 1 kHz are achieved with the ZlMPOL's specially designed CCD and are also quite high 
in the other instruments mentioned above. 

A FLC plate was chosen mainly for two reasons: low cost and the ability to switch 
fast. We wanted to explore the possibility of improving the accuracy of spectropolarimetric 
measurements using fast modulation. Pockels cells and piezoelastic modulators, both capable 
of switching at a rate of many kHz, were also examined but we found these devices to be 
too costly for our very limited budget. A piezoelastic modulator does not work at low 
frequencies and therefore would require a special masked detector which we do not have 
access to. Mechanical rotation would not allow us to explore the potential advantages of fast 
modulation and is also more expensive than a FLC plate. Another type of liquid crystal, 
liquid crystal variable retarders or LCVRs, could have been employed in a fast switching 
polarimeter. Such devic es are used, for examp le, in the HiVlS spectropolarimeter on the 
3.67-m AEOS telescope (JHarrington et al.ll2010l ). It should be noted that LCVR plates are 
almost two orders of magnitude slower than FLC plates. After considering several different 
options we decided that an FLC plate was the most cost effective option allowing us to 
explore fast modulation. 



2.3. Modulation 



The role of the half-wave plate is to switch the orthogonal states of linear polarization 
produced by the preceding quarter-wave plate. When the optical axis of the half-wave plate 
is aligned with the polarization plane the linearly polarized beam passes through the half- 
wave plate unchanged; when the optical axis of the half- wave plate is at an angle of 45° to the 
polarization plane it rotates the linearly polarized light 90°. As a result, the half- wave plate 
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effectively switches the two orthogonal polarizations by rotating both polarizations 90°. 

Because of the beam displacer that follows the half- wave plate (see below), switching the 
polarization states also results in the polarimeter's two output beams exchanging places. In 
dimaPol the modulation of the beams is synchronized with the CCD detector. As discussed 
above, the switching needs to be relatively fast in order to eliminate instrumental effects that 
can affect accurate polarization measurements. Simply reading out the detector frequently 
after short exposures is not desirable since this would significantly increase overheads because 
of the CCD readout time and likely also quickly wear out the instrument's shutter. The SITe 
detector is, however, a three-phase CCD that allows bi-directional charge transfer. This 
permits rapid shuffling of accumulated charge on the CCD back and forth between rows 
many times before the detector is read out. Such parallel transfers are almost three orders 
of magnitude faster than the read out time for the detector. Rather than reading out the 
CCD frequently we instead shuffle the charge on the CCD up and down with every switch of 
the half-wave plate. This charge shuffling technique for measuring polarization signals was 



first implemented by iMcLean et al.l (119811 ) in the Imaging SpectroPolarimeter at the Royal 
Observatory, Edinburgh. The shuffle distance in dimaPol is chosen to be the same as the 
final separation of the two output spectra, or about 40 pixels in a direction perpendicular to 
the spectra. During the charge shuffling the orthogonal polarization states never cross each 
other's path. The CCD shutter remains open for the duration of an exposure for multiple 
switch cycles and is kept closed during CCD readout. 

The switching time of the PLC half- wave plate at room temperature is about 70 /xs, 
and while it increases to about 120 /xs at +10 C, switching at kHz rates is possible even at 
lower temperatures. This is significantly faster than the approximately 2 ms that is needed 
to move the charge 40 rows on the CCD so there is some 'dead time' while charge is being 
shuffled. Since the wave plate switching time is only a small fraction of the shuffle time the 
switching is done in the middle of the charge shuffling to minimize the crosstalk between the 
two polarizations. 



2.4. Retardance and the switching angle of the half-wave plate 

An obvious advantage of using an electrically switchable wave plate is that there are no 
moving parts in the polarimeter. This provides better stability, reliability, and significantly 
reduces both construction and maintenance costs. PLC wave plates are also true zero order 
retarders and they have a large field of view. The acceptance angle of our PLC wave plate 
is 20°. This permits its use in the //18 beam from the telescope without introducing any 
noticeable error in wave plate retardance. One disadvantage of PLC wave plates is the fact 
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Fig. 3. — Measured retardance of the FLC half wave plate used in dimaPol. The dashed line 
is the plate's nominal A/2 retardance. 

that they are chromatic; retardance of the DAO FLC wave plate differs from A/2 by 0.015 A 
in the middle of the spectropolarimeter's spectral range and increases to 0.027 A and 0.022 A 
at the blue and red ends respectively (Figure |3]). 

The data shown in Figure H] and Figure E] were obtained in our optical laboratory. The 
plate to be tested was placed first between two open (the axes parallel) and then two closed 
(the axes orthogonal) BVO Vis22 linear polarizers from Bolder Vision Optik. The quarter- 
wave plate was oriented at 45° with respect to the first fixed polarizer. The change in the 
beam intensity was measured with an Ocean Optics USB2000 fiber optic spectrometer and 
the wave plate retardance was calculated using the equation 

'III -L 



(5 (waves) 



27r 



■ arccos 



III + In 



where /|| and /+ are the beam intensities with the polarizers open and closed respectively. 
The retardance of the FLC plate was measured in the same way. All of the measurements 
were carried out at room temperature. 

While the FLC plate retardance is not temperature dependent the switching angle is. 
To examine this effect we carried out an experiment in the cold chamber at the HIA where 
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the temperature of the wave plate could be controlled. The wave plate was placed at 45° 
between two open polarizers. The second polarizer was rotated from its "open" position until 
a minimum in the beam intensity was reached. The deviation angle of the second polarizer 
from the "open" position is then 2x the departure of the wave plate switching angle from 45°. 
One of our spare plates was used in this test and the measurements were performed at 4600 A 
where the retardance of the plate is closest to ideal. Figure H] shows the results. We find that 
the switching angle is exactly 45° at +23 ± 1 C, a typical summer temperature at the DAO. 
In winter Victoria temperatures drop to +5 C and the switching angle at this temperature 
is approximately 48.5°. Our modelling shows that this change in the switching angle can 
introduce crosstalk from linear to circular polarization of up to 13%. The crosstalk should be 
significantly less in summer. Crosstalk from linear to circular polarization is not uncommon 
in other polarimeters but does not pose a problem when measuring circular polarization in 
upper main sequence magnetic stars because their intrinsic linear polarization is at least an 
order of magnitude lower than their circular polarization. Thermal stabilization of the FLC 
could solve the crosstalk issue but has not been implemented for dimaPol because of lack 
of space and the limited impact such crosstalk has on the science we anticipate performing 
with the instrument. 



2.5. Crosstalk 

Departures of wave plate parameters from the ideal creates crosstalk between the two 
polarimeter channels. Some of the light from the ordinary beam leaks into the extraordinary 
beam and vice versa. As a result one will see a change in the measured Stokes V. We used 
the data in Figures [2] through H] to estimate how much the Stokes V measurement is affected 
by such crosstalk. We obtained estimates for two different temperatures: +5 C and +15 C. 
We assumed that the instrument is fed with 100% circularly polarized light (i.e., the Stokes 
V is equal to 1) and the results are shown in Figure [51 The departure of the Stokes V from 
1 tells how much crosstalk there is in the system. As expected, the amount of crosstalk is 
higher at the lower temperature mostly due to the increase in the FLC switching angle from 
47° at +15 C to 48.5° at +5 C. This crosstalk can be significantly reduced by rotating the 
FLC plate by 2° around the beam and in the direction opposite to the wave plate axis tilt. 
The amount of crosstalk is reduced by a factor of two at the lower temperature and decreases 
insignificantly at the higher temperature. As a consequence, in the case of a simple Zeeman 
triplet and a magnetic field aligned parallel to the line of sight dimaPol will underestimate 
the longitudinal field component by about 0.5 to 1%. However, longitudinal magnetic fields 
are very rarely measured with an accuracy su fficient to detect such a small discrepancy (for 



example, iBychkov. Bychkova. fc Madejl 119771 ) 
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Fig. 4. — The FLC switching angle as a function of temperature. The data points are 
measurements performed in a cold chamber at the HIA for one of our spare wave plates. 
The nominal switching angle of the wave plate is shown by the dashed line. The solid line 
is a fit to the data of the form ci + C2a/(To — T), where T is the plate temperature and ci, 
C2, and To ar e constants. Such a tempe rature dependence is typical of that observed in FLC 
wave plates (JBlinov fc Chigrinovlll994j ). 



The difference in transmittance of the FLC wave plate along the fast and slow axes 
(up to a few percent), together with different transmittance for ordinary and extraordinary 
beams inside of the spectrograph, results in spurious continuum polarization. The amplitude 
of this signal reaches 1% at some wavelengths but since this unwanted signal is much broader 
than any spectral line in the region, including H/3, and is a smooth function of wavelength 
it can be easily filtered out during data reduction as discussed below. 



2.6. Beam displacer 

A 20 mm long calcite beam displacer separates the orthogonal linear polarizations into 
two parallel beams. The ordinary beam with the plane of polarization along the dispersion 
travels straight through the displacer while the extraordinary beam emerges from the beam 
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Fig. 5. — Depolarization due to crosstalk between the two channels of dimaPol (open circles). 
100% input circular polarization is assumed (Stokes V equal to 1). The crosstalk can be 
significantly reduced by rotating the FLC plate by 2° (filled circles). 
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Fig. 6. — A schematic representation of one switch cycle of dimaPol. From left to right: 
expose, move the accumulated charge up, switch the FLC wave plate and expose again, 
move the charge back down and switch the wave plate in preparation for the start of the 
next switch cycle. As a result, three spectra form on the CCD (see text). 



displacer with a constant displacement of about 2 mm in the direction perpendicular to the 
dispersion. Because of this, after passing through the spectrograph the two beams arrive 
at two different locations on the detector where both orthogonal states of polarization can 
be measured simultaneously. At the detector the spectra are separated by ~ 0.6 mm or 
K, 40 pixels. Again, it is important to note that the two beams never move on the detector 
and that there are no moving parts in the polarimetric module. The two beams always 
illuminate the same rows of the detector. 

With no half- wave plate switching and no CCD charge shuffling the observation of a star 
with dimaPolwovXd obviously result in two spectra being imaged on the CCD, but with very 
poor polarimetric precision because of instrumental effects mentioned above. The implemen- 
tation of charge shuffling results in three spectra being produced as follows (and illustrated 
in Figure E]). At the start of an exposure the left circularly polarized light traveling through 
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the polarimetric module becomes the ordinary beam while the right circularly polarized light 
becomes the extraordinary beam. As a result two spectra, LO and RE, appear on the CCD. 
Then, when the FLC wave plate switches (and ignoring for now the fact that the charge 
shuffling takes considerably longer than the wave plate switch) the charge in the LO and RE 
spectra is shifted upward so that RE is brought up to where the ordinary beam is projected 
on the CCD and LO is placed 40 rows higher, outside of the light collecting area. But at the 
same time, because of the wave plate switch the right circularly polarized light becomes the 
ordinary beam (RO) and so it is projected right onto RE. RO and RE then mix together to 
form a single right circularly polarized spectrum. The left circularly polarized light becomes 
the extraordinary beam after the wave plate switch and because of the charge shuffle this 
creates a third spectrum (LE). Now there are three spectra on the CCD: LO, RO+RE, and 
LE. When the modulator switches again all three spectra are shuffled 40 rows back down 
and the cycle repeats. How these three spectra are processed and manipulated to measure 
a stellar magnetic field is discussed in the next section. 

Since only a single beam displacer is used in dimaPol, the ordinary and extraordinary 
beams travel slightly different optical path lengths and therefore have different foci. ZEMAX 
models show that the difference in the best focus position is very small, approximately 20 /xm. 
When the spectrograph is focused at an intermediate position (when both images are slightly 
out of focus) the degradation in image quality amounts to less than 2 //m or roughly 0.1 pixel. 
This is less than our typical focus uncertainty and can, therefore, be neglected. The difference 
in optical path between the two beams is also wavelength dependent and as a result the size 
of the extraordinary image along the dispersion is about 15% larger at the blue end of the 
detector. There is also a small difference in dispersion between the two beams. Despite these 
effects, when the ordinary and extraordinary beams are combined (see the next section) the 
spectral resolution is reduced by less than 5% and so we felt there was little need to use a much 
more costly Savart plate instead of a single calcite crystal. Savart plates are comprised of 
two plates of birefringent material rotated 90° with respect to each other. The extraordinary 
beam in the first plate then becomes an ordinary beam in the second plate and vice versa 
and, as a result, both beams have nearly identical optical paths and focus positions. 

After the beam displacer a fast mechanical Uniblitz shutter is installed. The shutter 
has opening and closing times of 18 ms and 12 ms respectively. 

All optical surfaces in the polarimetric module are antireflection coated. The transmis- 
sion of the polarimeter module is therefore quite high and varies smoothly between 74 and 
80% within its operational wavelength range. 
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3. Fast Switching and Observing Strategy 

We have aheady noted that switching a spectropolarimeter quickly helps to reduce 
the impact of spurious signals that can arise from changes in the spectrograph aperture 
illumination due to seeing and transparency variations, guiding errors, and other telescope 
instabilities. Some of these changes in illumination can be averaged out during long exposures 
but trends that manifest themselves only on time scales of hours are not compensated for 
in this way. It therefore seems that a better observing strategy for effectively eliminating 
the effects of changes in aperture illumination is to take multiple short exposures, switch 
the polarimeter module after each, and then combine the resulting sequence of exposures. 
When a short switching time is used it is much more efficient to do the averaging on-chip by 
moving the accumulated charge back and forth rather than after many CCD readouts and 
this is what we do with dimaPol. 

But the question arises as to how many times the polarimeter module needs to be 
switched in order to provide the most effective averaging and hence elimination of spurious 
signals. We have performed a series of on-sky experiments in an attempt to answer this 
question. In these tests we kept the total integration time on the target constant but varied 
the shuffle time or the number of shuffles per exposure. We define the shuffie time as the 
amount of time used to expose briefiy, switch the polarimetric FLC half-wave plate, and 
move the charge on the detector. The number of shuffies is simply the total exposure time 
divided by the shuffle time. Only an even numbers of shuffles are used to make sure that 
the polarized spectra obtained for successive targets always appear in the same place on the 
CCD. 

A cycle of two shuffies therefore consists of the following sequence: expose, switch the 
wave plate, move the charge, expose, switch the wave plate to its original state and move the 
charge back. As noted earlier, it takes approximately 2 ms to move the charge 40 rows on 
the CCD while the switching of the wave plate is considerably faster and is performed at the 
middle of the shuffiing phase. The total dead time per shuffie is, therefore, approximately 
2 ms and this sets a lower limit to the possible shuffie time to ensure an efficient duty cycle. 
We refer to this as 'dead time' since during this period some of the light on the detector 
gets spread in between the spectra and this signal is removed, together with the rest of the 
background signal, during processing (described below). 

For our switching tests we observed Vega which, along with its brightness, offers the 
advantage of having narrow spectral lines and hence makes it an easy target to obtain high 
quality magnetic ffeld measurements in a very short period of time. The exposure time for 
individual observations in a given observing sequence was kept constant at 6 s to provide 
approximately the same signal per exposure, while the shuffie time from observation to 
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observation was varied between 0.06 and 3 s (corresponding to a range of 2 to 100 shuffles 
per exposure). CCD charge tr apping sites begin to manifest them selves when a few hundred 



shuffles or more are executed (JBland-Hawthorn &: Bartonlll994j ). These traps create spikes 



in the extracted spectra, polluting the polarization signal. In order to avoid this effect 
we limited ourselves to 100 shuffles. On three different nights we obtained a total of 21 
magnetic field measurements derived from the H/3 line. Each measurement was based on 
the average of 30 individual 6 s exposures. We did two things to minimize the effect of 
changing weather conditions and flux variations on our results: fflst, the shuffle time was 
changed quasi-randomly from measurement to measurement but within the 0.06 to 3 s range, 
and second the longitudinal field error bars were normalized by the flux as a{Bi)noTm = 

a{Bi) y^ flux/ flux max- 

Figure [7] shows how the accuracy of magnetic field measurements obtained with dimaPol 
vary as a function of the number of shuffles. The magnetic error bar is largest when only two 
wave plate switches and charge shuffles are carried out. The gain in accuracy has therefore 
been measured by normalizing the estimated error in the longitudinal magnetic field for an 
arbitrary even number of shuffles to the mean value found for observations with just two 
shuffles. The magnetic field error bars get smaller as the number of shuffles increases until 
more or less a plateau is reached at about 30 shuffles. The accuracy is improved by a factor 
of 2 when 30 or more switches are executed. We also carried out observations of scattered 
solar light. The telescope was pointed at the sky near the zenith during the day and series 
of 50 short 4 s exposures were obtained for several different shuffle speeds. Unlike the stellar 
measurements, these scattered solar light observations are not affected by slit losses and can 
therefore be used to test for systematics within the spectropolarimeter itself, while excluding 
other telescope instabilities. The results are presented in Figure [7] along with the results for 
Vega. 95% confidence intervals of the standard deviation are also shown for both Vega and 
scattered solar light. For Vega a third degree polynomial approximation of points was used 
for A^ < 30 and the average accuracy gain of 2.08 was assumed for A^ > 30. The points are 
within the confidence intervals for both data sets. There is no change in the accuracy gain 
seen for the scattered solar light data which means that most likely slit losses are responsible 
for the accuracy degradation in Vega at low switching speeds. Switching more than 30 times 
eliminates the effect these losses have on the longitudinal magnetic field measurements. 

As a result of these tests, we have adopted a standard observing strategy for use with 
dimaPol. To make a single magnetic field measurement we obtain a series of typically 10 to 
16 short exposures (each, for example, 5 or 10 min long for a sixth or seventh magnitude star) 
of the target with 60 wave plate switches and charge shuffles performed during each exposure. 
This keeps the switching time shorter than the typical time of instrumental response changes, 
and also allows an estimate of the magnetic field uncertainty to be made from the scatter of 



o 

o 
o 

< 



-15- 



2.5 



1.5 



"1 I I I I L_L_|__I__I__I- 

/ 
/ 

/ • 



-l__l__l__|__l_ J_ J__| _ 



/ 

° /° ° O 

/...°....°.o 

L*J I I I I I I \ I I L 



J I I I I I L 







20 40 60 80 

Number of shuffles 



100 



Fig. 7. — Measured variations in the accuracy of longitudinal magnetic field measurements of 
Vega (filled circles) and scattered solar light (open circles) as a function of the number of wave 
plate switches and charge shuffles executed by dimaPol. The gain in accuracy is measured in 
comparison to the average estimated magnetic field error when only two shuffles are carried 
out. 95% confidence intervals of the standard deviation are also shown for Vega (dashed 
lines) and for scattered solar light (dotted lines). 



measurements obtained from the individual exposures. 

Because exposure times are short, we also bin every two pixels of the SITe CCD along 
the direction perpendicular to the dispersion to speed up the readout of the CCD and reduce 
the impact of readout noise. Since typical seeing at the DAO is around 2 — 3 arcsec, we use 
the 2 arcsec entrance aperture for all observations with the polarimeter. We also acquire 60 s 
iron-argon comparison lamp exposures before the start and at the conclusion of the observing 
sequence for each star. The polarization module switches at 1 Hz for the arc exposures so 
that we have the usual 60 switches and shuffles for the arc observations that we use for our 
science exposures. 
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4. Data Reduction 

We have developed a collection of MIDAS scripts that enable near-real-time processing 
of the data produced with dimaPol. This software provides a measurement of the longitudinal 
magnetic field and its associated uncertainty within a few minutes of the acquisition of the 
last spectrum in an observing sequence for a star. For every individual exposure we perform 
bias subtraction, cosmic ray and background removal, and then extract the spectra. The bias 
level in each image is determined by measuring just the CCD overscan area and subtracting 
this from all of the other image pixels. The remaining background is removed in the next 
step. 

As has been mentioned earlier, the CCD shutter remains open for many wave plate 
switch and charge shuffle cycles. Because of this, even though the parallel charge transfer 
takes very little time, some light from the source does get spread out between the three 
recorded spectra during the shuffling process and this creates an additional background 
source on the detector. This extraneous light needs to be removed before the polarized 
spectra can be extracted. The background flux has a trapezoid like shape when viewed 
along the CCD columns. We divide each CCD image into three sections, one for each of the 
three spectra, and treat the background in each section separately. In each of the two outer 
sections we fit second order polynomials to the background flux measured on either side of 
each spectrum; a simple flrst order polynomial is used to average the background flux on 
either side of the spectrum in the middle section. These three independent polynomials are 
flt to every CCD column and the resulting background flt for each column is then subtracted 
from the corresponding column in the original image. 

The next processing step is the extraction of the polarized spectra. The three spectra 
(denoted LO, RO+RE, and LE above) are extracted separately. During the spectrograph, 
CCD and polarimeter conflguration carried out in advance of each polarimeter run, consid- 
erable care is taken to ensure that the acquired spectra run parallel to the CCD rows. The 
spectrum inclination and curvature are then typically less than 5 /xm over the approximately 
26 mm length of CCD. This makes the spectral extraction process quite simple: CCD rows 
inside of a predeflned aperture are simply co-added. We have determined that the best re- 
sults are obtained when flve CCD rows are coadded. Since the CCD is binned by a factor of 
two perpendicular to the dispersion this corresponds to ten unbinned rows on the detector. 
Such an aperture gives the best signal-to-noise ratio when a 2arcsec spectrograph circular 
entrance is used. The size of the aperture is flxed and the same size is used to extract all 
three spectra. The three apertures are placed 40 pixels from each other (corresponding ex- 
actly to the charge shuffle distance) to make sure that we always extract the same pixels for 
the RO-I-RE, LO, and LE spectra. The two polarizations are exposed on exactly the same set 
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of physical pixels even though the charge during shuffling may be stored in different places. 
After extraction the LO and LE spectra are combined into a single spectrum, LO+LE. 

As alluded to earlier, the PLC wave plate and other optical elements in the spectrograph 
introduce an extremely broadband spurious continuum polarization that can produce a weak 
circular polarization signal, even from an unpolarized light source. The effect is wavelength 
dependent and for an unpolarized light source leads to an LO+LE spectrum that is ~ 
1% brighter than the RO+RE spectrum at both ends while the two spectra have almost 
identical intensities in the middle. Even though this broadband signal changes very little from 
exposure to exposure and from object to object, it can affect precise line shift measurements 
required for magnetic field observations since the small wavelength-dependent difference in 
brightness between the two polarizations results in an artificial line shift. It is also dependent 
on the size of the spectral window used for measuring pixel shifts and discussed below: the 
broader the spectral window the larger the introduced pixel shift. If not taken into account 
during the data processing, this instrumental shift can be as large as a few hundredths of a 
pixel which, as will be seen, can lead to a spurious, and obviously significant, signal on the 
order of a few 100 G. 

In order to remove this effect we apply a chromatic correction to our data. A third 
degree polynomial is sufficient to fit the variation in /lo+le/-^ro+re along the direction 
of the dispersion. The unnormalized RO+RE spectrum is then multiplied by this fit to 
eliminate the artificial brightness difference between the two polarizations, and removes the 
spurious line shift as well. Note that the chromatic correction is only derived from and 
applied to the final, combined LO+LE and RO+RE spectra for an entire observation and 
not to individual exposures. This minimizes the possible impact of low signal-to-noise levels 
in individual observations on the polynomial fit for the correction. We feel that normalizing 
-^lo+le/-^ro+re is far preferable than trying to normalize the intensities of the LO+LE 
and RO+RE spectra independently. The spurious continuum polarization as it manifests 
itself in /lo+le/-^ro+re can be fit very well with a low order polynomial. On the other 
hand the continuum in the intensity spectrum is considerably more complex in shape and is 
more difficult to reconstruct, especially when broad spectral lines are present and when the 
normalization has to be carried out in a completely automated fashion for stars with a wide 
range of spectral types during data acquisition. 

Finally, the magnetic shift observed in a single spectral line due to the Zeeman effect is 
then measured by performing a Fourier cross-correlation of the final combined and corrected 
LO+LE and RO+RE spectra in a spectral window centered on the spectral line of interest. 
Several different methods of measuring the magnetic shift in a single line have been used by 



different authors. These include the center of gravity method (JRees &: Semel Ill979[ ) and its 
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impleme ntations bvlMathva (11988 



and iDonati et al. j 1997 ), the so-cal 



ed ph otopolarimetric 



Borra &: Landstreetl (1973) for use with a 



Bagnulo et all (l2002[ ) 



method (iBray fc Loughheadl Il964l ) adapted by 

two-channel Balmer-hne photoelectric Pockels cell polarimeter and by 

for use with the FORSl spectrop olarimeter, and the Fourier cross-correlation technique 

( jMonin. Fabrika. fc Valyavinll2002l ). We have adopted the Fourier cross-correlation method 

for use with dimaPol because of the technique's relative insensitivity to continuum and line 

shapes which permits its use with unnormalized spectra and without prior knowledge of the 

intrinsic profile of the line in the absence of a magnetic field. 

The magnetic shift is translated into a longitudinal magnetic field, Bi (G), a ccording to 
the li ne's magnetic sensitivity or the Lande factor and the relation (for example, iLandstreet 
1992h 

Bi = /\X X 0.15/(2 X 4.67 x lO^^^A^c/), 

where AX is the observed shift between the spectra in pixels, A is the wavelength of the 
lin e center (A) and q is the Lande facto r for the line. For hydrogen lines the Lande factor is 
1 (jCasini fc Landi Degl'Innocentilll994 ) while for metallic lines the Lande factors are taken 
from the VALD database (JKupka et al.ll2000l ). The measured shift in pixels is multiplied 
by the average dispersion of 0.15 Apixel"^. The change in dispersion that arises along 
the spectrum as a result of the calcite beam displacer is small and can be neglected. For 
H/3 observations with dimaPol this equation is transformed to the simple relation Bi = 
6.8 X AXkG. In other words, a one pixel shift between the RO+RE and LO+LE spectra 
corresponds to a longitudinal magnetic field strength of 6.8 kG. 

More than ten individual observations, and hence magnetic measurements, are usually 
obtained in succession for a typical V = 7 target and many more are often acquired for ob- 
servations of very bright stars (usually known magnetic 'standards' or 'nulls' observed during 
each observing run to monitor performance of the instrument). The multiple independent 
longitudinal field measurements are then averaged and an error bar can then be derived from 
the scatter of these individual values, after weighting by the flux as discussed in the previous 
section. 

Figure M shows a histogram of 100 measurements of the shift between the RO+RE 
and LO-I-LE spectra for the bright, sharp-lined Ap star HD 112413 (a^ CVn) obtained in 
succession over almost one hour on a single night. The mean shift between spectra is ~ 
0.09 pixels corresponding to a field strength of ~ 600 G and the sample standard deviation is 
~ 0.04 pixels or ^ 270 G. The measurements are apparently distributed normally as a curve 
representing a normal distribution with a mean of 0.09 pixels and variance of (0.04 pixels)^ 
fits the histogram of observations very well. A Kolmogorov-Smirnov test confirms that the 
measurements are very likely normally distributed: the P- Value of 0.98 indicates that the 
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Fig. 8. — An illustration of the apparent normal distribution of individual magnetic field 
observations obtained with dimaPol. The histogram shows the distribution of 100 individ- 
ual measured shifts (in pixels) between left- and right-circularly polarized spectra obtained 
during a single observing sequence of the bright, magnetic Ap star a^ CVn. The superim- 
posed curve shows an expected normal distribution with a mean of 0.09 pixels and variance 
of (0.04 pixels)^. A Kolmogorov-Smirnov test confirms that the observed pixel shifts are 
indistinguishable from a normal distribution. 



histogram points and the normal distribution differ insignificantly. When the points are 
normally distributed the standard error of the mean is proportional to iV~^/^ where N is 
the number of observations. Therefore, it seems quite reasonable in the case of dimaPol 
data to estimate the magnetic field error bars by multiplying the standard deviation of the 
individual observations by this same factor. In the above example the error bar is then 
0.04 X 100-1/2 ^ 0.004 pixels or ^ 27 G. 

Experimentation has demonstrated that the error bars are minimized by optimizing the 
size of the window used for the Fourier cross-correlation so that the window includes the 
part of the profile that gives the strongest polarization signal (Figure [H]). The window size 
ranges from 16 pixels or about 2.4 A for narrow-line stars (such as a^ CVn) to 128 pixels or 
19 A for broad-line stars. 
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Fig. 9. — An example of dimaPol observations of the H/3 line of the prototypical magnetic 
Ap star a^ CVn. The Stokes V signature is at the bottom and has been multiplied by a 
factor of 10 relative to the / spectrum. Uncertainties in longitudinal field measurements 
derived from such data are minimized by optimizing the size of the spectral window used 
for the Fourier cross-correlation used to measure the shift between the LO+LE and RO+RE 
spectra by ensuring that it only includes the portion of the line profile that gives the strongest 
polarization signal. For a^ CVn the spectral window is indicated by the vertical dashed lines 
and is 16 pixels wide. For more rapidly rotating stars the window can be up to 128 pixels 
wide, or about 19 A. 

A careful reader may think that we have omitted some details of our data reduction 
procedure but this is not the case. Flat fielding of the DAO spectropolarimetry is not 
necessary for our primary goal of measuring longitudinal magnetic fields. As mentioned 
earlier, both of the left- and right-circularly polarized spectra are exposed on exactly the 
same CCD pixels. As a result, pixel-to-pixel sensitivity variations are the same for both and 
do not affect the final magnetic line shift determination. We do, however, normally obtain 
sequences of high signal-to-noise fiat field spectra several times during each observing run 
in case they are needed in the future, for example, for a detailed analysis of Stokes I line 
profile variations. 



The use of a Fourier cross-correlation algorithm to derive the observed longitudinal 
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magnetic field also simplifies processing by removing any need for continuum normalization. 
Thanks to the fast switching of the polarimeter module, variations in spectrum brightness 
due to changing sky conditions, slit losses, and other instrument changes are effectively 
averaged for both polarizations. This results in the shape of the continuum being almost 
identical in both polarizations. The Fourier cross-correlation algorithm that we use is insen- 
sitive to features that have a similar shape in both polarizations and so the algorithm can 
be successfully applied to unnormalized spectra which greatly facilitates the near-real-time 
processing of the polarimeter data while observing. 

Wavelength calibration of the spectra is also not required for our magnetic field mea- 
surements since the Fourier cross-correlation of the polarized spectra is carried out in pixel 
space. Wavelengths are generally only needed for line identification and for this we use the 
two iron-argon comparison lamp exposures obtained before the start and at the conclusion 
of each sequence of observations for a given science target. The comparison spectra are 
processed in the same way as the science spectra, the arc lines are identified, and a second 
order polynomial is used to derive the dispersion relation. 

In addition to providing measurements of longitudinal magnetic fields, dimaPol also 
allows investigations of the circular polarization signal within the profiles of many spectral 
lines, including H/3. The circular polarization signal, Stokes V , is calculated very simply as 
the difference in brightness between the LO+LE and RO+RE spectra, normalized by the 
total flux: 

y = (-^LO+LE — -^RO+Re)/(-^LO+LE + -^RO+Re)- 

An example of a short section of a Stokes V spectrum can be seen in Figure [91 

The presence of the spurious continuum polarization introduced by the FLC wave plate 
and other optics and the crosstalk from linear to circular polarization makes it very difficult 
to measure circular continuum polarization with dimaPol. However, in sources with a high 
degree of circular polarization and an insignificant amount of continuum linear polarization 
this can be done, albeit with somewhat limited precision. 



5. Magnetic Field Observations with dimaPol 

Details of observations of the first new magnetic star discovered with dimaPol have 



recently been presented by iBohlender fc MoninI (1201 ll ). Along with such ongoing science 



programs, however, regular testing of the polarimeter is also performed to confirm that the 
instrument does not produce spurious polarization signals that could lead to false magnetic 
field detections. To ensure that this is the case, during every scheduled polarimeter run we 
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observe magnetic 'null' standards, bright stars that have been observed with high precision 
numerous times with a variety of other polarimeters and have shown no evidence for de- 
tectable magnetic fields at a level of a few tens of G. Seven such stars have been repeatedly 
observed by us over a two year period. These include HD886 (7 Peg), HD 36486 (5 0ri A), 
HD 76644 (iUMa), HD 97633 (^Leo), HD 156164 (5 Her), HD 172167 (Vega), and HD 210027 
(iPeg). 

Both slowly and rapidly rotating stars are included in this collection of standards, with 
t^sini ranging from 3kms~^ to 290kms~^. 133 individual magnetic field measurements 
have been obtained from observations of the H/3 line and they are all within ±3cr of OG. 
The average error bar, asi is 78 G and the best is 21 G. The null standards do not show 
any systematic errors at this level and data for some of them, shown in Figure dUl also 
demonstrates that the measured longitudinal magnetic field error bars scale with the square 
root of the flux as expected. Using the single H/3 line, an accuracy of 20 to 40 G is reached 
in 20min in good weather conditions for our slow rotating (< 20km/s) third to fourth 
magnitude standards. Obviously the accuracy is lower for faster rotating stars simply because 
of the decrease in sensitivity of the cross-correlation measurements as the spectral lines 
become broader. For example, for HD 76644 {V = 3.1; vsini = 150km/s) we obtain 70— 80 G 
error bars in about 30 min. For the typical sixth to seventh magnitude stars in our science 
programs an accuracy of 100 to 200 G is reached in between one and two hours. 

Vega has recently been shown to possess an extremely weak magnetic field (< 1 G, 



Lignieres et al.ll2009l : iPetit et al.ll2010l ) but a field of this level is well beyond the detection 
capability of dimaPol and therefore we continue to use this convenient, bright star as a 
magnetic null standard. As shown in (Figure [TTj) . no evidence for a magnetic field is seen in 
the Stokes V signature of Vega as observed with dimaPol. 

One null standard, iPeg, has been observed 21 times on 15 different nights. The scatter 
of the resulting 21 individual H/3 longitudinal magnetic field measurements about the G 
level is 48 G, in good agreement with the average ctbi = 58 G for the individual observations 
(Figure [T2|) . This again reassures us that our derived error bars adequately reflect the true 
accuracy of the observations. 

As an additional test of dimaPol performance we have carried out observations of scat- 
tered solar light during several observing runs. For such observations the telescope is usually 
pointed at the sky near the zenith during the day and long series of short exposures (typi- 
cally 4 to 6 s, although as long as 30 s on several occasions) are obtained. Very high precision 
measurements can be acquired because of the brightness of the light source and the high 
contrast provided by the solar spectrum's narrow lines. We have carried out these observa- 
tions 37 times over a period of two years and a total of 3660 individual spectra have been 
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obtained. No systematic shifts between the individual RO+RE and LO+LE spectra have 
been detected at a level corresponding to a longitudinal magnetic field strength of 15 G. 

Our discussion above shows that dimaPol does not produce spurious magnetic field 
detections at least at the level of a few tens of G. This does not, however, guarantee that the 
instrument can reliably measure a real magnetic signal when a field is present. In order to 
regularly test the instrument's ability to detect and accurately evaluate longitudinal magnetic 
fields we also attempt to make observations of a few well-known magnetic Ap stars during 
each observing run. One example of an object which we have repeatedly observed is the 
bright, prototypical magnet ic Ap star HP 112 413 (a^ CVn). The star has a well-established 
rotation period of 5^46939 (JFarnsworthl 119321) and a strong, reversi ng longitudinal magnetic 
field that varies between -1 and -|-lkG (JKochukhov fc Wadd 120101 ). At the time of writing, 
we have observed it 34 times with total integration times ranging from 7min to one hour, 
with an average of 20 min. A representative observation showing the magnetic signature in 
the Stokes V profile is illustrated in Figure [131 Note that the scale used for the Stokes V 
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Fig. 10. — The uncertainty of the longitudinal magnetic field as measured with dimaPol as 
a function of source fiux. Data for the Sun and two magnetic null standard stars are shown: 
stars - scattered solar light, filled points - iPeg, open points - 7 Peg. The solid lines show 
the yj flux fit for each source. 
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spectrum in this figure is the same as that used in Figure [TT] for Vega; the signature of the 
magnetic field in o? CVn is obvious. 

Longitudinal magnetic field measurements of o? CVn have been obtained by using both 
the H/3 line and the single metallic line Fe II A4923. The latter feature was selected because 
it is the strongest metallic line in the wavelength range covered by dimaPol and, with a 
Lande value of 1.69, it is also quite sensitive to the presence of a magnetic field. Table [1] 
provides the details of the magnetic field observations including the HJD at the mid-point of 
the spectropolarimeter data acquisition, the longitudinal magnetic field values derived from 
both spectral lines, and the rotation phase derived from the ephemeris of lFarnsworthl ( 119321 ): 
JD = 2419869.720 + 5.46939 x E. 

The resulting magnetic field curves produced from the measurements of the H/3 and 
Fe II lines are shown in Figures [T41 and [T^ respectively. Figu re [T^ compares dimaPol H;9 lon- 
gitudinal magnetic field measurements to those published by lBorra fc Landstreetl (119771 ) and 
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Fig. 11. — A typical dimaPol observation of the magnetic null standard star Vega. The 
Stokes / spectrum is shown at the top of the figure while the lower spectrum is the Stokes 
V signature multiplied by a factor of 10 relative to the / spectrum. The total exposure time 
of the observation was 3 min and the longitudinal magnetic field strength derived from the 
H/3 hue is 8 ± 22 G. 
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Fig. 12. — dimaPol magnetic field measurements obtained from H/3 observations of the 
magnetic null standard t Peg over an approximately two year period. The scatter of the in- 
dividual measurements around G is in good agreement with the mean size of the individual 
error bars, confirming that the error bars are reasonable estimates of the uncertainties in 
the observations. Note that three additional data points have been omitted from the plot; 
all three were obtained on a cloudy night and have large error bars (within 3o" of G) that 
make inclusion in the plot difficult. 



obtained with a two-channel Balmer-line photoelectric Pockels cell polarimeter . The Pockels 



cell p olarimeter measurements have been multiplied by 4/5 as suggested by iMathys et al. 



(I2OOOI ) in order to account for the Stark effect. This correction only applies to measurements 
in the line wings where the Stark effect is a factor. The agreement between the two sets of 
observations is very good. Note that the error bars for our new DAQ data a r e, on average, 
almost two times better (56 G versus 98 G) than those of iBorra fc Landstreetl (119771 ) despite 
the approximately three times shorter integration time for the new observations. The scatter 
of the DAO points around a best fit sinusoid is 60 G, quite close to the average error bar of 
our measurements. This is another indication that our error bars adequately represent the 
accuracy of our magnetic field measurements. 

Figure [15] shows our longitudinal magnetic field measurements obtained with the Fe II 
A4923 line as well as values derived from the Least Squares Deconvolution (LSD) of multiple 
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metallic lines (jWade et al.ll2000[ ) using spectra obtained with the MuSiCoS spectropolarime- 
ter at Pic du Midi. Again, the two sets of measurements are in very good agreement. We 
also again point out the fact that the average size of the error bars of our single-line magnetic 
field measurements (40 G) is slightly better than that of the multi-line LSD measurements 
(55 G) obtained with a comparable telescope. 

As noted earlier, one of the science goals for which we designed dimaPol was an investi- 
gation of the occasional differences observed in the longitudinal magnetic field measured for 
a single star with the H/3 line versus metallic lines. Comparison of Figures [T^ and [13] clearly 
illustrates such a difference for magnetic field measurements of a^ CVn and a detailed dis- 
cussion of this (occasionally very large) effect in this and other stars, and its possible origin, 
will be presented in a future paper. 
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Fig. 13. — A typical dimaPol observation of the prototypical magnetic Ap star a^ CVn. The 
Stokes I spectrum is shown at the top of the figure while the lower spectrum is the Stokes 
V signature multiplied by a factor of 10 relative to the / spectrum. The total exposure time 
of the observation was 15 min and the longitudinal magnetic field strength derived from the 
H/3 line is -965 ± 41 G. 
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Fig. 14. — A comparison of longitudinal magnetic field observations of a^ CVn obtained 
from H/3 line observations and phased on its 5?46939 rotation period. The solid points are 
new o bservations obtained with dimaPo/ while the open points are from lBorra fc Landstreet 
(119771 ) after correction for the Stark effect (see text). The solid curve is a best-fit sinusoid 
to the DAO data. 



6. Conclusion 

We have designed and constructed a very inexpensive dual-beam polarimeter for use 
on the DAO 1.8-m Plaskett Telescope. The polarimeter, dimaPol, is implemented as a 
plug-in module for the telescope's Cassegrain spectrograph and provides medium resolution 
circular spectropolarimetry of point sources over a 260 A wide spectral window centered on 
the H/3 line. The instrument therefore provides DAO 1.8-m observers the ability to measure 
longitudinal magnetic fields using the Zeeman effect in both the H/3 and metal lines of the 
peculiar A and B stars of the upper main sequence. 

We demonstrate that the use of an electrically switchable ferro-electric liquid crystal 
(FLC) half-wave plate to quickly switch the orthogonal polarizations between the instru- 
ment's two channels, in combination with synchronized shuffling of the accumulated charge 
on the CCD, greatly reduces instrumental effects and increases the accuracy of the magnetic 
field measurements. With extensive on-sky tests we show that by executing 30 to 100 such 
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Fig. 15. — As in Figure [H] but for metallic line observations. The solid points are new 
observations obtained with dimaPol from the single Fe II A492 3 line while the op en points 
are measurements derived from LSD of multiple metallic lines ( jWade et al.ll2000l ) observed 
with MuSiCoS. The best sinusoidal fit to the DAO H/3 measurements (the same as in Fig- 
ure \\M is also shown. Note a significant difference between the H/3 fit and the metallic line 
measurements near the negative extremum (see text). 



switches the accuracy of the longitudinal magnetic field measurements can be improved by a 
factor of two. The use of the FLC wave plate also eliminates the need for any moving parts 
in the polarimeter module and therefore provides better stability, reliability and significantly 
lower construction and maintenance costs. 

A near-real-time data reduction pipeline has been implemented with MIDAS and pro- 
vides observers with the ability to obtain a longitudinal magnetic field measurement and its 
uncertainty within seconds of completing an observing sequence on a star. The observed shift 
between the recorded left- and right-circularly polarized spectra is performed with a Fourier 
cross-correlation technique that is not sensitive to continuum variations or intrinsic spectral 
line shapes and so can be carried out without interaction from the user. For magnetic field 
measurements using the H/3 line dimaPol provides the simple relation Bi = 6.8 x AXkG 
where AX is the observed shift between the polarized spectra in pixels. 
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A number of bright, well-established magnetic stars as well as magnetic 'null' standards 
have been repeatedly observed over a two-year period. No systematic errors have been found 
at the level of 15 to 20 G. Our multiple observations of the magnetic Ap standard o? CVn (as 
well as other objects not discussed here) clearly demonstrate that the longitudinal magnetic 
field derived from dimaPol H/3 and metal line observations agree well with previous published 
results. In addition, the derived uncertainties of the field measurements compare favorably 
with or exceed the accuracies achieved with other often much more expensive polarimeters 
on similar sized telescopes. To provide a rough estimate of the capabilities of dimaPol, 
observations of a typical sixth or seventh magnitude star in our current science programs 
reach an uncertainty of 1 00 to 200 G in a typical one to two hour observation (e.g., see 



Bohlender fc Moninll201lh . 



We have begun several research programs with dimaPol on the DAO 1.8-m telescope 
but encourage any interested users to apply for observing time with the instrument. Our 
main survey program consists of an extensive search for magnetic fields in previously poorly 
studied (and generally relatively faint, V > 7) Ap and Bp stars. In a second program we are 
taking advantage of the capability of the instrument to provide longitudinal magnetic field 
measurements in both the H/3 line as well as metal lines to investigate the reality and cause 
of the often very discrepant magnetic field strengths measured with the two diagnostics, but 
frequently with very different instrumentation. As discussed in the previous section, a"^ CVn 
is an example of one such object. 

Despite its cost of only several thousand dollars, dimaPol is proving to be a very ca- 
pable addition to the complement of instrumentation available on the DAO 1.8-m Plaskett 
Telescope and has opened up an entirely new avenue of research at the DAO. At such a low 
cost many institutes with spectroscopic capabihties on even very modest telescopes could 
readily consider building a similar device for research purposes or student training. 

We would like to thank J.D. Landstreet, G.A. Wade, G.A. Chuntonov, and V.G. Shtol' 
for very valuable discussions. J.H. Grunhut also provided valuable help during the first 
dimaPol commissioning runs. 

Facilities: DAO:1.85m. 
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Table 1. Longitudinal magnetic field measurements for the Ap star c^ CVn obtained with 

dimaPol. 



HJD 

(-245 0000) 


Bi (G) 

(H/3) 


Bi (G) 
(Fell) 


4> 


4480.98881 


-394 ± 


41 


-342 ± 


19 


0.1771 


4487.04114 


185 ± 


54 


183 ± 


35 


0.2837 


4487.98377 


685 ± 


45 


711 ± 


26 


0.4561 


4543.89436 


-95 ± 


54 


-140 ± 


34 


0.6785 


4632.77958 


-965 ± 


41 


-694 ± 


62 


0.9299 


4633.74205 


-704 ± 


34 


-596 ± 


27 


0.1059 


4634.71760 


282 ± 


61 


213 ± 


57 


0.2843 


4634.73870 


337 ± 


69 


211 ± 


47 


0.2881 


4636.78101 


-14 ± 


50 


21 ± 


27 


0.6615 


4640.74014 


681 ± 


76 


539 ± 


34 


0.3854 


4659.76615 


-841 ± 


39 


-706 ± 


23 


0.8640 


4660.77011 


-817 ± 


41 


-641 ± 


35 


0.0476 


4661.74017 


-41 ± 


50 


-195 ± 


39 


0.2250 


4870.00537 


311 ± 


40 


236 ± 


22 


0.3033 


4877.95372 


-421 ± 


64 


-531 ± 


47 


0.7565 


4878.08744 


-561 ± 


48 


-627 ± 


34 


0.7810 


4878.85740 


-937 ± 


60 


-778 ± 


42 


0.9217 


4879.03109 


-1013 ± 


33 


-708 ± 


42 


0.9535 


4902.93131 


307 ± 


50 


454 ± 


30 


0.3233 


4903.93935 


623 ± 


72 


651 ± 


40 


0.5076 


4934.75652 


-555 ± 


47 


-507 ± 


34 


0.1421 


4936.94704 


618 ± 


28 


614 ± 


16 


0.5426 


4937.92873 


-149 ± 


50 


-379 ± 


32 


0.7221 


4939.84380 


-785 ± 


44 


-616 ± 


27 


0.0722 


4939.86512 


-802 ± 


59 


-591 ± 


39 


0.0761 


4959.77708 


-201 ± 


68 


-287 ± 


45 


0.7168 


4960.88101 


-930 ± 


37 


-666 ± 


30 


0.9186 


4963.73155 


684 ± 


99 


668 ± 


67 


0.4398 


4964.76458 


157 ± 


63 


209 ± 


52 


0.6287 


5058.69226 


-836 ± 141 


-628 ± 


84 


0.8020 


5059.69697 


-986 ± 


65 


-746 ± 


50 


0.9857 


5068.69036 


160 ± 100 


269 ± 


61 


0.6300 


5321.81108 


-965 ± 


47 


-662 ± 


50 


0.9095 


5322.78558 


-701 ± 


51 


-576 ± 


39 


0.0877 



